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Abstract 

Electromagnetically induced transparency-like behavior in the extreme ultraviolet (XUV) is stud- 
ied theoretically, including the effect of intense 800 nm laser dressing of He 2s2p (^P°) (^P°) and 2p^ 
(^S'^) (^S'') autoionizing states. We present an ab initio solution of the time-dependent Schrodinger 
equation (TDSE) in an LS-coupling configuration interaction basis set. The method enables a rig- 
orous treatment of optical field ionization of these coupled autoionizing states into the N = 2 
continuum in addition to = 1. Our calculated transient absorption spectra show encouraging 
agreement with experiment. 
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5 I. INTRODUCTION 



The use of two laser fields to dress the eigenstates of a three-level system has attracted 
significant attention over past two decades The interaction of light with a coherent 
atomic ensemble leads to important quantum interference phenomena, including electro- 
magnetically induced transparency (EIT) , in which the laser-induced coherence of 
atomic states leads to quantum interference between the excitation pathways that control 
the optical response. It can be used to manipulate the linear and nonlinear susceptibilities 
of matter from the near-infrared to ultraviolet. This has resulted in applications such as 
slowed and stored light as well as nonlinear frequency conversion {4]. 

Presently there is great interest in the application of light in the X-ray regime, produced 
by high-order harmonics and free electron lasers, to investigate novel coherent X-ray optical 
phenomena. In this context, ultrafast electron dynamics of the autoionization processes 
in atoms has gained attention in recent years. The advent of femtosecond and attosecond 
XUV laser pulses in a pump-probe setting has made it possible to follow these processes 
in real time js]. These experiments involving two coherent laser pulses allow study of the 
time-resolved interaction of the direct ionization pathway with pathways involving one or 



21 more resonant states. 



his opens the possibility of new ways to monitor and control the au- 



toionization dynamics p, l7| and consequently the XUV absorption properties of the dressed 
medium. 

n n 

Theoretical studies due to Bachau et al. [8[ and Themelis et al. [9] predicted the transfor- 
mation of the original unperturbed Fano line-shape 10|] into an Autler-Townes doublet due 
to Rabi flopping between a pair of coherently-coupled doubly-excited states. Loh et al. 
observed EIT-like behavior in an XUV probe induced by coherent coupling of the 2s2p (^P°) 
and 2p^ (^S°) doubly-excited states in helium, where the XUV probe was created from IR 
laser-produced high-order harmonics. Specifically, that study used a femtosecond infrared 
(IR) pulse both to couple the doubly-excited states and to produce the probe pulse. Buth 



et al. 12| predicted an EIT for X-rays in laser dressed neon. More recently, Ranitovic et al. 



13| used a two-color multiphoton ionization of helium in combined IR laser and XUV high 



harmonic fields to control the transparency. They showed that this scheme can induce full 
electromagnetic transparency. 

Quantum mechanical calculations of the photoabsorption spectra for the dressed atoms 



36 represent a challenge for theory, especially when the dressing laser pulse couples doubly- 

37 excited states. The time-consuming solution of the time-dependent Schrodinger equation 

38 (TDSE) is required for every photon energy of the probe pulse. An alternative theoretical 

39 treatment has been introduced by Tong and Toshima [14] and successfully used in Ref. [3]. 

40 It is based on the Floquet representation and on the calculation of the autocorrelation 

41 function that is indepdendent of the probe laser photon energy. Therefore, the time and 

42 memory demanding numerical calculation of the autocorrelation function is performed only 

43 once. Since the photoabsorption cross section has a simple analytical dependence on the 

44 autocorrelation function, its calculation is not difficult. However, the approach introduced in 

45 Ref. fy| involves a diagonalization of the Hamiltonian matrix in the Floquet representation, 

46 typically having dimensions 10^ x 10^ and this step is still challenging for a supercomputer. 

47 To our knowledge, there is no method which enables a routine theoretical calculation of 

48 the photoabsorption spectra in the presence of a strong dressing laser pulse. The aim of the 

49 present work is to make a step towards the improvement of the theoretical treatment of this 

50 problem. This study deals with the modified Fano line-shape in the photoabsorption spec- 

51 trum of the He 2s2p (^P°) doubly-excited state caused by the coupling of this state with the 



52 2p^ (^S*^) doubly-excited state by an IR femtosecond pulse. Loh et al. jll| studied this effect 



53 experimentally and interpreted obtained results as EIT-like behavior. A simplified three- 



54 level [15| model treatment of this problem was included in Ref. The model describes 

55 formation of an Autler-Townes doublet and takes into account the effect of the optical field 

56 ionization of the doubly-excited states due to the IR dressing laser using a somehow ad hoc 

57 adaptation of the Ammosov-Delone-Krainov (ADK) method. Our approach presented in 

58 this work is based on the numerical solution of the TDSE. It takes into account the pos- 

59 sible coupling of more than two autoionizing states and the continuum in principle more 



accurately than the approach presented in Ref. Uj. Therefore, the present treatment of 



61 the optical field ionization of the doubly-exited states due to the coupling laser field is more 

62 quantitative as well. Photoabsorptio n sp ectra calculated using parameters corresponding to 

63 the experimental setup used in Ref. ll| are compared with results obtained for reduced IR 

64 pulse intensity and for significantly reduced IR detuning. Comparison of all these results 

65 allows us to study the effect of the optical field ionization and of the IR laser detuning on 

66 the formation of EIT-like structures in the spectra. 

67 The rest of this paper is organized as follows: The theoretical approach and assumptions 



68 of our model are discussed in Sec. HIl A detailed description of the field free Hamiltonian 

69 and static photoabsorption spectrum is given in Sec. IIIIl while Sec. [IV] discusses the pho- 

70 toabsorption yields of the dressed atom. 

71 Atomic units are used throughout the rest of the paper. 



II. THEORETICAL APPROACH 



73 The mode 

74 Loh et al. 



introduced in the present study is based on the experimental publication of to 
Uj. In that experiment the 2s2p (ip°) (ip°) and 2p2 (^S^) (^S^) doubly-excited 

75 states of He are coupled by a strong 800 nm field. The 2s2p (^P°) and 2p^ (^S*^) states 

76 are located at 60.15 [l^ and 62.06 eV 17] above the Is^ ground state, respectively. This 

77 system is then probed after a variable time delay by an XUV pulse. The XUV probe pulse 

78 is produced with a table-top, laser-based setup by high-order harmonic generation (HHG). 

79 The authors use a commercial Ti:sapphire laser system (2.4 W, 800 nm, 42 fs, IkHz) to 

80 produce the optical dressing and HHG beam. The estimated photon flux at the source is 10^ 
31 photons per pulse for the high-order harmonic centered at 60.2 eV. The XUV pulse duration 
82 is estimated to be 30 fs FWHM. The dressing pulse peak intensity is 1.4 x lO^^W/cm^ 



83 and its duration is 42 fs FWHM. The spectrometer used in [ll| measures the transmission 

84 spectrum instead of the absorption spectrum. The unknown fraction of dressed atoms in the 

85 interaction volume prevents the photoabsorption spectrum of the dressed atom from being 

86 retrieved experimentally. Therefore, authors measure the transient absorption spectra and 

87 identify EIT-like structures. 

88 The XUV laser pulse couples the Is^ ground state with the 2s2p (^P°) autoionizing state 

89 with energy 35.56 eV above the = 1 ionization threshold. Therefore, it is reasonable 

90 to assume that every absorption of an XUV photon leads to photoionization. The aim of 

91 the present work is to calculate the population of the ground state at the end of the pulse 

92 as a function of the XUV photon energy hwx for different dressing-probe time delays to- 

93 The computed loss of the ground state population is interpreted as the total photoabsorp- 

94 tion yield. Its comparison with the yields obtained when the IR pulse is absent allows for 

95 recognition of EIT-like structures. 

96 Helium is a two-electron system. The two-particle basis set used in the present calcula- 



97 tions consists of L5'-coupled independent-particle basis functions 



18 



19| . The Hamiltonian 



98 for the interaction of each electron with the core is 

Mr) = -iv^ - p (1) 

99 where r is the coordinate of the electron with respect to the nucleus. The corresponding (ra- 

100 dially rescaled) single-particle radial wave function (f{r) for an electron with orbital angular 

101 momentum I, is given by 

1 0^2 /(/ + !) 2\ . , , , , , 

+ (2) 



2 (ir^ 2r^ r\ 

102 The atomic Hamiltonian including the interaction between the electrons is 

// = Mri) + Mr2) + — , (3) 

103 where ri2 = |ri — is the distance between the electrons. 

104 To calculate the time development of the ground state population we solve the TDSE for 

105 the two linearly-polarized laser fields in the length gauge: 

i^*(ri, r2, t)^[H + Di(t + to) + Dx(t) + Vc] *(ri, ra, t), (4) 

106 where Di^x{t) — E/_x(^) • d = E/_x(^) • + i"2) are the interaction operators for the IR 

107 and XUV laser pulse in the dipole approximation, respectively. Vq denotes the complex 

108 absorbing potential (CAP) added to treat ionization properly within the context of a finite 

109 volume calculation. As in the experimental convention, a positive time delay to means that 

110 the XUV probe pulse center arrives later than the pump pulse center. The electric fields 

111 ^i{t) and Ex(i) are given by 

Ei,x{t) = zEi,x{t) = zElxfit) cos {uji,xt) , (5) 

112 where cuj^x are the IR and XUV laser field frequencies, z is the unit vector along the po- 
ns larization axis, E^ ^ are the peak amplitudes of the electric fields and fi,xiS) are the pulse 
114 envelopes. Throughout this work it is assumed that both envelopes have the cos^ form 



hit) = I 



cos^ (yrt/r/) , -Ti/2<t <Ti/2 

VolScos^ (7r(i - tp)/Tip) , -Tip/2 <t-tp< tip/2 (6) 
0, otherwise. 



115 resp. 



cos^ {Tit/rx) , -Tx/2 < t < Tx/2 
Jx[t) = \ {7) 
0, otherwise, 

116 where r/ x denote the durations of the respective laser pulses. The second term in Eq. ([6]) 



117 denotes the IR postpulse which appears in the experiment 



111 ] with time delay t 



118 and duration tjp = 30 fs. The CAP used in the present calculations has the form 20 | 



88 fs 



0, ri < n and r2 < n 

Vc — \ (oj 
—irj [(ri — Th)'^ + (r2 — rf,)^] , otherwise, 

119 where ?7 > is the strength of the CAP and is the inner boundary of the absorber. The 

120 TDSE (jlj) is integrated in a box using the expansion of the time-dependent wave function 

121 \E'(ri, r2, t) in a basis set consisting of eigenstates of the field-free atomic Hamiltonian H: 

imax N'l 

^(ri,r2,t) = J2Y.CLnmLniri,r2), (9) 

L=0 n=l 

122 where i"2) denotes the n**^ eigenfunction of H having total angular momentum L and 

123 corresponding to the real "box" eigenenergy em'- 

H^Ln{ru ra) = eL„$L„(ri, r2). (10) 

124 The expansion ([9]) includes Nl eigenstates for the angular momentum L, and it includes all 

125 the angular momenta up to Lmax- The basis set used in the variational calculations of the 

126 field-free eigenstates <l'L„(ri, consists of LS-coupled independent-particle basis functions. 

127 Antisymmetric two-electron basis functions coupled to form a state of definite S,L, and 



128 parity P are constructed in terms of one-electron orbitals ([2]) as explained in Refs. [18|, Il9 |. 

129 The radial one-electron wave functions ifji{r) are obtained by solving Eq. ([2]) in a radial box 

130 of size ro > with boundary conditions (Pji{0) = iPji{ro) = for both positive and negative 

131 energies e. The index j denotes the number of nodes of the radial wave function inside the 



2l| 



132 box. Eq. ([2]) is solved by expansion of the radial orbitals in the basis set of B-splines 

133 Bf{r) satisfying the bound-state boundary conditions: 

i,k 

134 The radial box is divided into several intervals (denoted by index k) by knots and the 

135 B-splines are defined in terms of the knots in each interval. 



136 Use of the LS-coupled basis as discussed above is nowadays a standard technique fre 



quently employed in the theoretical treatment of both time- dependent [19|, |21|, |22| and 
time-independent 18[ problems for two active electron atomic systems. The radial one- 



19 



21 



22| or discrete variable 



electron orbitals are often expressed in terms of B-splines 
representation (DVR) basis sets |23 |. 

In order to suppress the reflections from the outer boundary of the radial box in the 
solution of the TDSE and to keep its size reasonably small, the CAP term (|H]) is added 
to the Hamiltonian used for the time-propagation. The presence of this term makes the 
total Hamiltonian of the system a complex symmetric rather than hermitian operator 20 1. 
One note should be made at this point regarding the form of the matrix elements of all the 
operators used here. Since all the calculations are performed in the eigenrepresentation of the 
field-free atomic Hamiltonian H not including the CAP, the operators H, Dj{t) and Dx{t) 
have real matrix elements, as it is convenient in hermitian quantum mechanics. Vc is the only 
non-hermitian term in the total Hamiltonian. Therefore, it is not necessary to distinguish the 
standard hermitian scalar product and complex symmetric product throughout the present 



work (see Ref. 



24| and references therein for details). 



Both laser pulses are linearly polarized along the 2;- axis, and of course the ^S*^ ground state 
of He has zero projection Ml of the total orbital angular momentum on the polarization 
axis. Accordingly it is also enough to consider only the 2;-component of the dipole operator. 
According to the selection rules for the matrix elements of the dipole moment operators 



{^Ln \ Zi + Z2 I^L'n')' thcsc are non-zero only if L — L' = ±1 [25]. 



Projection of the TDSE (jl]) on the eigenstates of the field-free Hamiltonian (UOl) yields 
the following set of coupled differential equations of the first order: 



ic{t) = [H + Vc + Dj(t + to) + T>x{t)] c(t). 



(12) 



159 where c(t) is the vector of expansion coefficients used in ([9]) and the matrices on the right 

160 hand side are the representations of the corresponding operators in the basis set $i„(ri, r2) 

161 introduced above. The dimension of this system of differential equations is = X^l^cT -^l- 

162 In the present representation H is a real diagonal matrix of the field-free atomic eigenen- 

163 ergies. Since Vc given by Eq. ^ has no angular dependence, Vc is an imaginary and 

164 symmetric block diagonal matrix with Lmax + 1 non-zero blocks having dimensions A^^:- 
166 D/ x(^) are real symmetric banded matrices with only non-zero blocks corresponding to the 



166 coupling of the total angular momentum L with L ± 1, because we only consider Ml = 

167 states in this study. The physical picture of the phenomenon studied here suggests that the 

168 IR dressing laser pulse couples the singlet autoionizing states 2s2p (^P°), 2p^ (^S°) and the 

169 continuum states made accessible by the IR pulse. The energy difference between the ground 

170 state and the lowest odd-parity, singlet excited state - ls2p(^P°) of the field- free atom is so 

171 high that it requires more than 18 photons to be absorbed to excite the ground state by 

172 the IR pulse. Transitions of higher orders would thus be necessary to excite higher excited 

173 states directly from the IR pulse acting on the ground state. Therefore, it is a reasonable 

174 assumption to neglect the IR coupling of the ground state with excited states. Note that the 



175 weak effect of the IR field for intensities considered in this work is also suggested in Ref. [26 1. 

176 Therefore, the matrix elements of the dipole operator in the interaction matrix for the IR 

177 field (d^)ooL'n' were explicitly set to zero (considering the ground state the first element in 

178 the basis set). This assumption allows for a dramatic improvement in the convergence of the 

179 solutions of the TDSE with respect to the size of the basis set. On the other hand, it is the 

180 primary role of the XUV probe laser pulse to couple the ground state with the autoionizing 

181 state 2s2p (^P°) and with nearby continuum states. Since it has rather high photon energy 

182 and low intensity (perturbative regime), it is unlikely that the XUV coupling between the 

183 doubly-excited states and the continuum states will play any significant role. Therefore, 

184 in the present study all the dipole matrix elements in the XUV interaction matrix Dx(t) 

185 except those coupling to the ground state (d^)ooL'n', are explicitly set to zero. Note that 
ar approximation for the XUV coupling in a similar context has been employed in 
i]. Using this approximation the system of coupled differential equations (fT2l) can be 

188 written as follows: 

tc{t) = [H + Vc + Ej{t + tD)di + Ex{t)df] c(t). (13) 

189 The structure of the matrices is showed in Fig. [H The atomic system is initially in the 

191 ground state. Considering the ground state to be the first element of the basis set, the 

192 initial condition can be written as Co,i(to) = 1; CL,n>i(i^o) = 0. The neglect of the ground 

193 state coupling with the other field-free eigenstates via the IR laser pulse allows for reduction 

194 of the total interval for the propagation of the TDSE as well. Since tj > tx, it would be 

195 necessary to start the propagation of the TDSE (fT2|) in the time to = ~tD — t//2, if this 

196 assumption were not made. It would also then be necessary to perform the propagation of 

8 
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FIG. 1. Structure of the matrices in the system of coupled equations (jl3p 
The grey areas denote the non-zero blocks corresponding to definite total angular momenta and 
their dipole coupling. Matrix dl, has zero first line and first column, while matrix has 

non-zero only the first row and column. 



the wave function until the time ti = max{— t/j + tj/2,tx/2}. Although the total wave 
function described by the vector c{t) changes at times t > ti due to the presence of the 
CAP, the population of the ground state will remain constant after the end of the pulses, 
since the ground state wave function is well localized in the region where the CAP is zero. 
Therefore, for the purpose of the present project it is enough to perform the propagation 
until ti and extract the population of the ground state co,i(ti). The approximation that 
the IR field does not affect the ground state allows for the start of the propagation at time 
to = —Tx/2 and terminate it at time ti = tx/2. In view of the long duration of both pulses, 
this is an important simplification of the numerical calculations. 

The method of solution of the TDSE used in the present work is based on the split 
operator method in the basis set representation and it is described in the appendix in detail. 
This represents an alternative approach to the higher-order quadrature methods 19l. 1271. |28| 



frequently used in context of the interaction of atoms and molecules with strong laser fields. 
Application of the split-operator technique is straightforward in case of the wave-packet 



211 propagation (for details see Ref. 29(| and references therein). However, the exponential of 



the time dependent operator remains a challenge in the basis set representation . Palacios 
et al. 30| introduced an implementation in the DVR basis set which has properties that 



allow for an efficient numerical calculation of these exponentials. In the implementation 
used here these time-dependent exponentials are calculated in the diagonal representation 
of the dipole operators, as is discussed in the Appendix. 

The photoabsorption yields obtained from the solutions of the TDSE ( ITSll together with 
the static photoabsorption yield enable a calculation of the transient absorption spectra. 
This allows to relate the theoretical results presented here with the experimental transient 



220 absorption spectra published in Ref. 



Hi 



221 III. FIELD-FREE HAMILTONIAN AND STATIC PHOTOABSORPTION SPEC- 

222 TRUM 

223 Although we have found that the eigenrepresentation of the atomic Hamiltonian H with- 

224 out the CAP is more practical for the solutions of the TDSE f[T^ . it is useful to study the 

225 spectrum of the complex symmetric Hamiltonian He = H + Vc- 

226 This representation allows for easier optimization of the CAP parameter r] and of the size 

227 of the radial box tq. The complex eigenenergies allow for a better insight into the structure 

228 of the discretized continuum and evaluate its quality for the solution of the full TDSE as 

229 well. The eigenfunctions 0i„(ri,r2) of the complex symmetric operator He = H + Vc are 

230 given by the equation 

Hc(j)Ln{ri, rs) = ELn(j)Ln{ri, Tg), (14) 



231 where the indices L,n have identical meaning as in Eq. ffTOj) . The complex eigenergies can 

232 be written as E^n = -S£„ — i^Ln/'^- They should reproduce well the energy positions i?£„ 

233 and widths F^^ of the 2s2p (^P°) and 2p^ (^S'^) autoionizing states for L = 0, 1. In addition, 

234 the density of the discretized continuum states must be sufficient to accurately describe the 

235 interaction between the autoionizing states and the continuum and to allow for converged 

236 calculation of the photoionization yield. These properties of the basis set are controlled by 

237 the size of the box tq, the position of the absorbing boundary rf, and the strength of the 



238 CAP ?7, as is described in Ref. 20j. These parameters have been optimized in the present 

239 calculations to the values tq = 200 a.u. and = 150 a.u. These values allow for a good 

240 representation of the complex eigenenergies of the autoionizing states, they yield negligible 

241 decay widths of the ground and low excited states, and they provide a sufficient number of 

242 continuum states above the ionization threshold. The radial box is spanned by 700 B-splines 

243 of the sixth order. The LS-coupled expansion includes the single-particle partial waves up to 

244 / = 19, every partial wave includes 200 radial orbitals (see Eq. ([2])). In order to describe the 

245 autoionizing states of the interest properly, the CI expansions of the states with L = and 

246 L = 1 include the doubly-excited configurations, where both electrons can occupy orbitals 

247 with / < 6 and j < 6. In the remaining configurations one electron is restricted to occupy 
243 a Is, 2s, or 2p orbital, while the other electron can occupy any orbital, constrained so that 

10 



249 the configurations contribute to the states with given spin 5", total angular momentum L 

250 and parity P. The higher partial waves I included in the expansion are required for the 

251 representation of the continuum states having high total angular momentum L; they are 

252 essential for obtaining converged solutions of the full TDSE including the nonperturbative 

253 IR laser pulse. Optimization of the parameter rj as is described in the Ref. \2Q^ using this basis 

254 set yields the value rj = 10~^. The complex energies of the autoionizing states calculated 

255 using the basis set presented here are compared with previously published experimental 

256 results 



IgI IitI in Table HI The widths of both autoionizing states calculated here are in good 



TABLE I. Comparison of the positions and widths of the 2s2p (^P°) and 2p^ (^S'') doubly-excited 
states calculated from diagonalization of the Hamiltonian He presented in this work with experi- 
mental values. The value r] = 10~^ was used for the CAP parameter. The positions are relative to 
the ground state Is^. The values in brackets are obtained by fitting the perturbation calculations 
to the Fano line-shape. 

Position (eV) Width (eV) 

Present calculation Experiment Present calculation Experiment 

2s2p (ip°) 59.687 (59.688) 60.15 [16] 0.0411 (0.0456) 0.0373 [16] 

2p2 (iS^) 61.747 62.06 \17} 0.0063 0.0059 [17] 

257 

258 

259 agreement with the experimental results. The discrepancies in their relative positions with 

260 respect to the ground state are mainly due to the inaccurate representation of the Is^ ground 

261 state by our truncated atomic Hamiltonian H. It gives a computed ground state energy of 



262 Eis-2 = —2.88343 a. u., while the accurate experimental value is Ef^^ = —2.90378 a. u. 31 ]. 

263 To improve the ground state energy (and subsequently the relative positions of the double 

264 excited states) a significantly larger basis set would be necessary. However, the difference of 

265 doubly-excited state energies is more important for the calculations of the photoabsorption 

266 spectra than their relative position with respect to the ground state. As can be seen from 

267 Table [J the energy difference according to the calculations is AE = 2.06 eV, while the 

268 experimental results give AE^^^ = 1.91 eV. 

269 The distribution of the ^S*^ and continuum states in the complex energy domain is 

270 plotted in Fig. [2j The IR laser pulse assumed here allows for the optical field ionization of 

11 
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FIG. 2. (Color online) Complex eigenenergies of the complex symmetric field-free Hamiltonian 
He spanning the interval of energies probed by the XUV laser pulse. Energies of states with the 
total angular momentum L = and L = 1 are plotted. The autoionizing states are labeled by 
their principal independent particle CI configurations. The autoionizing state with energy 62.22 eV 
without label does not have a single dominant CI configuration, its discrete component is a linear 
combination of the 2s3p and 3s2p configurations. The N = 2 ionization threshold is at energy 
64.86 eV. 

I — I 

271 the doubly-excited states into the N = 2 continuum. Since Loh et al. [11] suggested that 

272 inclusion of this effect is important for the proper treatment of the phenomenon studied 

273 here, it is important to have a sufficiently high density of continuum states with higher total 
278 angular momenta L in the energy range above the N = 2 ionization threshold. Present 

276 calculations of the dipole matrix yields fiab = 2.11 a. u. for the matrix element coupling 

277 the 2s2p (^P°) and 2p^ (^S'^) states. This is in a good agreement with previously published 



278 value 2.17 a.u. obtained from the eigenchannel R-matrix calculations 

279 To understand the effect of the IR laser pulse it is a necessary prerequisite to calculate 

280 the static photoabsorption spectrum in the absence of the dressing laser pulse. The static 

281 photoabsorption yield corresponds to a Fano line-shape [lOj across the known energy ui of 



282 the 2s2p (^P°) doubly-excited state with respect to the Is^ ground state, which has a known 

12 



FIG. 3. (Color online) Loss of the ground state population in the absence of the IR pulse and its 

comparison with the Fano line-shape convolved with the bandwidth of the pulse. 
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width Fi and value qi of the Fano hne-shape parameter. However, it is a stringent test of 
our basis set quahty and of the CAP used throughout this work to calculate this photoab- 
sorption yield independently using the present Hamiltonian and compare the corresponding 
spectrum with previously published results ll|. This can be done by setting = in 
Eq. f|T3l) . Taking into account the low intensity of the XUV laser field, it is most tractable 
to calculate the static photoabsorption spectrum using time-dependent perturbation theory. 
The corresponding loss of population of the ground state is shown in Fig. [3l This curve 
should correspond to the convolution of the Fano line-shape with the the XUV laser pulse 
bandwidth. Therefore, we have fitted the absorption yield obtained from the perturbation 
calculations to this convolution and extracted the position, width and line-shape parameter 
of the 2s2p (^P°) autoionizing state. The fitting procedure assumes a gaussian envelope of 
the laser pulse with bandwidth 0.06 eV (FWHM). The values for position and width obtained 
from the fitting are compared with previously published results in Table [H As expected, 
they agree very well with the values obtained from diagonalization of He- In addition, the 
fitting procedure yields qi = —2.71, in good correspondence with the value -2.75 published 
in Ref. ll|]. The convolved Fano line-shape corresponding to the fitted parameters is also 



shown in Fig. [3l The good correspondence between the perturbation calculations (using the 
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FIG. 4. (Color online) Loss of the ground state population for several time delays as a function 
of the XUV photon energy calculated for the IR laser pulse with Air = 800 nm and peak intensity 
= 1.4 X 10^3 W/cm^ 

301 cos^ XUV pulse envelope) and the Fano line-shape fit convolved with the gaussian function 

302 suggests that the results are not sensitive to the details of the XUV pulse envelope. 



303 IV. PHOTOABSORPTION BY THE DRESSED HELIUM ATOM 

304 The loss of the ground state population at the end of the probe pulse p{uJx) = 1 — |coo(''^i) |^ 

305 calculated from Eq. ( IT3l) as a function of the XUV photon energy ux for several time delays 
308 t£) is shown in Fig. HI Since the duration of the dressing laser pulse is longer than the 

308 duration of the XUV pulse, different intensities of the IR field are probed at different time 

309 delays to. As can be seen in Fig. HI the Fano line-shape disappears in the case of the perfect 

310 overlap of both pulses and the photoabsorption yield becomes essentially structureless. The 

311 peak and the minimum shifted with respect to the static case become more pronounced 

312 with increasing time delay t^- An additional peak appears around the XUV photon energy 

313 60.6 eV for > Ofs. Its energy position does not change significantly with changing time 

314 delay. Loh et al. performed a model calculation of the photoabsorption yield based on 
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315 the three-level model developed in Ref. 15| describing the Autler-Townes doublet formation. 



316 It explains the presence of the additional peak in the photoabsorption yield as well as the 



317 shift of the higher peak with respect to the static case. The Ref. model calculations 

318 that include optical field ionization of the doubly-excited states into the N = 2 continuum 

319 due to the IR field confirm that this effect leads to a significant line-width broadening. This 

320 broadening makes the Autler-Townes doublet nearly invisible and the position of the peak 

321 at positive XUV detunings hardly shifts as the IR pulse intensity is changed. Since the 

322 calculations presented in this work were performed at an intensity that allows for optical 

323 field ionization, our results are in good correspondence with the theoretical considerations 



324 presented in Ref. 



325 The present calculations account for the time profile of both pulses, whereby the pho- 

326 toabsorption yields shown in Fig. H] together with the static photoabsorption yield plotted in 

327 Fig. [3l enable a calculation of the transient absorption spectra that can be compared with 



328 experimental results published in 



329 The change of the optical density AOD is calculated using the formula that takes into 



330 account the energy resolution of the spectrometer used in the experiment HI] : 



I\UU{ujx) - - logio 7 r — -„ — ^273"' \^^) 

j exp [-ao{u)nti - [u-ujx) \duj 

331 where (t^u) and aoioj) are the photoabsorption cross section of the dressed medium and 

332 static photoabsorption cross section in the absence of the IR pulse, respectively, rtt is the 

333 atom number density of the helium gas in the tube and ^ is the length of the tube in 

334 the direction parallel to the direction of laser propagation. Auj = 0.18 eV is the energy 

335 resolution of the spectrometer (FWHM) used in the experiment pjj]. Determination of the 

336 photoabsorption cross section from the photoionization yield obtained from Eq. (fT3|) is not 

337 straightforward. The usual relation between the single-photon absorption cross section and 

338 the photoabsorption rate (probability per photon per second) dp/dt for an infinitely long 

339 monochromatic pulse with photon energy u is given by the formula 

= ^ J, (16) 

340 where dp/dt is the photoabsorption rate as a function of the photon energy u and / is the 

341 corresponding light intensity. This formula assumes that the photoabsorption rate does not 

342 depend on time. The situation in the present work is more complicated. Both pulses have 
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343 a finite duration, and the photoabsorption rate varies with time. This time dependence 

344 should be included in Eq. f lT6|) . However, when one takes into account other assumptions 

345 made regarding the shape of the XUV pulse, it is a good approximation to calculate the 

346 photoabsorption cross section using the formula 

a{ux) = A—— — , (17) 

347 where p{ujx) is the photoabsorption probability obtained from Eq. f[T51) . Tx is the duration 

348 of the XUV pulse and Ix is the corresponding XUV pulse intensity used to calculate p{ujx)- 

349 The scaling factor A approximately accounts for the time dependence of the XUV field 

350 intensity and of the photoabsorption rate. Its value was chosen so that the calculated 

351 absorption spectra reproduce well the experimental transient absorption spectra published 



352 in Ref. Since A is related to the XUV pulse only, the value used to calculate the 



353 static optical density is the same as the value used to calculate the optical density in the 



354 presence of the IR field. The comparison between our theoretical results and the Ref. [11 1 

356 experiment is presented in Fig. [51 The best correspondence between the present calclations 

357 and experimental results was achieved by setting A = 3.5 for the time delays to = Ofs, 20 fs 

358 and 40 fs. The remaining spectra with time delays to = QOk and 80 fs compare best with the 

359 experimental results for A = 5. Therefore, the best fit of the experimental data is obtained 

360 when the product IxTx used in the calculations is smaller than the values estimated in the 

361 experimental work [11]. The calculated transient absorption spectra plotted in Fig. [5] have 

362 been shifted on the energy axis by 0.45 eV to correct for the inaccuracy of the 2s2p (^P°) 

363 doubly-excited state position with respect to the ground state in the present calculations (see 

364 Table [T] and the discussion in Sec. IIIII) . A change of the value of the parameter A affects the 

365 ratio of the peak magnitudes in the transient absorption spectra, however it does not affect 

366 the positions of the structures in the curves. As can be seen in Fig. [5l the present theoretical 

367 calculations are generally in encouraging agreement with the experimental results. The left 

368 peak that develops with increasing time delay is a consequence of the rising maximum of 

369 the photoabsorption spectrum due to decreasing IR intensities probed by the XUV pulse. 

370 The figure shows that this peak develops in the present calculations more rapidly with 

371 increasing time delay than is evident in the experimental results. The peak at positive XUV 

372 detunings calculated in the present study corresponds well to the experimental results. Since 

373 the probe pulse is produced by the HHG, it is expectable that its spatial and frequency 
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ated using Eq. (fT5]l for different time 



FIG. 5. (Color online) Transient absorption spectra calcu 
delays t£) and their comparison with experimental results 1]|. The value of the scaling constant 
A = 3.5 was used to calculate the spectra with time delays Ofs, 20 fs and 40 fs. The value A = 5 



was used to calculate the spectra with time delays 60 fs and 80 fs. 



374 profile have more difficult structure than a simple gaussian envelope as it is assumed in this 



375 work. Therefore, t 

376 published in Ref. 



le discrepancies between the experimental transient absorption spectra 



11| and the spectra presented here can possibly be attributed to the 



377 simplified shape of the XUV pulse assumed in the present theoretical treatment. 

378 Although the photoabsorption yields in Fig. H] show the reduced XUV photon absorption 

379 caused by the dressing IR laser pulse, the relation between these results and the EIT-like 

380 structures is not straightforward. The spectral peaks corresponding to the dressed states 

381 of the Autler-Townes doublet approach each other linearly as the Rabi splitting Qc of the 

382 states coupled by the IR laser pulse decreases. This behavior is not apparent in Fig. SI 

383 The Rabi splittings Qc probed by the peak of the XUV laser pulse at time delays 1^ = fs, 

384 20 fs and 40 fs are flc = 1.18 eV, 0.62 eV and 0.0025 eV, respectively. The splitting for larger 

385 time delays tu is significantly smaller than the width of the 2s2p (^P°) doubly-excited state. 

386 However, Fig. H] does not show any corresponding rapid change in the energy separation of 
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FIG. 6. (Color online) Photoabsorption yield calculated for the peak intensity of the IR field 
/q = 7 X 10^^ W/cm^ and for the dressing pulse wavelength A/ = 600 nm corresponding to the 
energy difference between the 2s2p (^P°) and 2p^ (^S'^) doubly-excited states. Curves for time 
delays t/j = Ofs, 20 fs 40 fs are compared with the static photoabsorption yield with no IR pulse 

387 the peaks as the time delay is changed. The clear interpretation of these photoabsorption 

388 yields is complicated by the optical field ionization of the doubly-excited states, by the 

389 finite durations of the laser pulses, and by the large detuning of the dressing laser from the 

390 energy difference of the 2s2p (^P°) and 2p^ (^S*^) doubly-excited states. In order to eliminate 

391 the effect of the large IR detuning, another calculation was performed with the wavelength 

392 of the IR laser pulse modified to 600 nm. The corresponding photon energy 2.06 eV now 

393 matches the calculated energy difference between the doubly-excited states 2s2p (^P°) and 

394 2p^ (^S'^) (see Table [T]). In order to suppress the effect of the optical field ionization, the peak 

395 intensity of the coupling IR laser pulse was decreased to 1° = 7x 10^^ W/cm^. Corresponding 
39S photoabsorption yields are shown in Fig. El This graph shows for to = Ofs two peaks with 

398 energy separation 0.2 eV and minimum between them at energy close to the position of the 

399 static photoabsorption peak. The Rabi splitting of the autoionizing states coupled by the 

400 IR laser pulse with the reduced peak intensity used here is Vtc = 0.27 eV. The curve for time 

401 delay to = 20 fs shows similar behavior. The energy separation of the peaks is 0.16 eV, the 
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402 Rabi splitting of the 2s2p (^P°) and 2p^ (^S*^) doubly-excited states probed by the peak of 

403 the XUV laser pulse is Qc = 0.145 eV. In case of larger time delay = 40 fs both peaks 

404 have a strong overlap and their energy separation is smaller than the width of the 2s2p (^P°) 

405 doubly-excited state. This correspondence confirms that we can interpret these structures 

406 as the Autler-Townes doublet. Energy separations of the peaks for single time delays to 

407 do not exactly match the corresponding Rabi splittings Qc, mainly because of the finite 

408 duration of both pulses. Another reason could be the coupling of additional doubly-excited 

409 states by the IR laser pulse that would add further complications to the three-level picture. 

410 The presence of an additional autoionizing state in the coupling scheme also appears to 

411 explain the smallest additional peak at an energy below 59.4 eV which occurs at small time 

412 delays in Fig. El The 2s^ (^S'^) doubly-excited state is a good candidate, since the energy 

413 difference between the 2s^ (^S*^) and 2s2p (^P°) doubly-excited states is 2.373 eV in the 

414 present calculations. In general, the calculations performed with the reduced IR intensity 

415 and with the photon energy adjusted as discussed above show the EIT-like structure more 

416 clearly than the calculations performed with parameters of the IR laser pulse used in the 

417 experiment. This suggests that the Autler-Townes structure in the calculations carried out 

418 using the experimental parameters are not very pronounced simply because the dressing 

419 laser detuning is so large. 

420 In order to investigate the manner in which the optical field ionization of the autoionizing 

421 states affects the dressed photoabsorption spectra, another calculation with the adjusted 

422 IR laser pulse (A/ = 600 nm) was performed. The experimental peak intensity of the IR 

423 dressing laser pulse Ij = 1.4 x 10^'^W/cm^ was used. The resulting photoabsorption yields 
428 for several time delays to are shown in Fig. [71 This graph shows a broad peak around 

426 the energy 58.7eV and a richer structure in the energy interval from 59.7eV to 60.3 eV 

427 in case of perfect temporal overlap of the two laser pulses. These structures approach 

428 each other with increasing time delay tr,. Since optical field ionization broadens the line- 

429 widths, the characterization of the individual peaks in the photoabsorption spectra is not 

430 straightforward. However, the well pronounced separation of the broad peak at negative 

431 XUV detuning and the structure at positive XUV detuning is most likely related to the 

432 coupling of the 2s2p (^P°) and 2p^ (^S'^) autoionizing states, as can be assumed from the 

433 change of the separation with increasing time delay tr,. It is possible that the IR field 

434 at the experimental peak intensity 7° = 1.4 x lO^^W/cm^ enables coupling of additional 

19 




FIG. 7. (Color online) The photoabsorption yield calculated for the peak IR field intensity Jo = 
1.4 X 10^^ W/cm^ and for the dressing pulse wavelength A/ = 600 nm corresponding to the energy 
difference between the 2s2p (""^P") and 2p^ (^S'^) doubly-excited states. Curves for time delays 
t£) = Ofs, 20 fs 40 fs are compared with the static photoabsorption yield. 

435 autoionizing states and this coupling leads to the rich structure of peaks appearing at positive 

436 XUV detunings. 



437 V. CONCLUSIONS 

438 This work presents theoretical treatment of the photoabsorption in XUV for helium 

439 dressed by the IR laser. The 2s2p (^P°) and 2p^ (^S*^) doubly-excited states are coupled by 

440 strong 800 nm laser pulse with duration 42 fs. This system is probed by the weak XUV pulse 

441 produced by HHG. The photoabsorption yields in the presence of both pulses are calculated 

442 as a function of the XUV photon energy for several dressing-probe time delays by solving the 

443 TDSE. The calculated transient absorption spectra are in reasonable correspondence with 

444 experimental results of Loh et al. The presence of the coupling IR laser significantly 

445 reduces the peak photoabsorption yield and changes of the time delay between the dressing 

446 and probe pulse show the dependence on the IR field intensity that is probed by the XUV 
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447 photons. Calculations performed with the wavelength of the dressing laser modified to 

448 A/ = 600 nm, corresponding to zero detuning of the IR laser, clearly show a structure in 

449 the photoabsorption spectrum that can be interpreted as the Autler-Townes doublet and an 

lat this was not clearly pronounced in the 



11| primarily because of the large detuning 



450 ElT-like structure. This allows us to conclude t. 

451 experimental results published in Ref. Loh et al. 

452 of the IR laser. The photoabsorption yields calculated for different dressing laser intensities 

453 showed that the optical field ionization of the doubly-excited states broadens the linewidths 

454 and increasing dressing laser intensity leads to formation of more complicated structures in 

455 the photoabsorption spectra. These could result from the increasing importance of additional 

456 autoionizing states at higher IR intensity. 

457 The method of solution of the TDSE presented in this work proved to be stable and 

458 suitable for theoretical treatment of two active electron atoms in pump-probe settings of the 

459 laser pulses with a potential of further improvement in the future. 
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467 Appendix: Method of the solution of the TDSE 

468 Since it is necessary to integrate the system of equations (fT2|) for a set of different XUV 

469 photon energies hujx and pump-probe time delays t£,, it is important to employ an integra- 

470 tion method that efficiently meets these requirements. After tests of several different meth- 

— I I 1 

471 ods used elsewhere in similar contexts (e.g. Run ge-K utta (l9| and Crank- Nicholson 29|) we 

472 chose to implement the split-operator technique |29| in a basis set representation. The basic 

473 idea of this propagation method is that if the total time-dependent Hamiltonian HT(t) can 

474 be written as a sum of the time- independent and time-dependent parts Hr(t) = Ho + Hi(t), 
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the propagation of the wave function between times t and t + At, where At is sufficiently 
small, can be implemented to second order in the time step, as follows: 



*(t + At) = exp 



If we denote 



iAt. 



exp 



t+At 



iii{t')dt' 



exp 



iAt^ 



Fi,x{t,At) 



t+At 



Ei,x{t')dt', 



Ho ) *(t) + O [{Aty] . 

(A.l) 

(A.2) 



then the application of this method to the present system of equations f lT^ yields 

iAt 



c(t + At) = exp 



c 



exp [-idiFi{t + to, At) - zdf Fx(t, At)] 
zAt 



exp 



H + V 



c) 



c(t) + O [(At)'] . (A.3) 



The first and third exponential terms in this equation can be calculated rather easily, as they 
are time independent (for a fixed time step At) and all the matrices are block diagonal. Since 
the eigenrepresentation of the field-free Hamiltonian H is used, this term can be expressed 
using the Crank-Nicholson formula 



exp 



iAt 



(H + V 



c) 



1 - i(H + Vc) At/4 



+ O [(At)^] , 



(A.4) 



1 + « (H + Vc) At/4 

which is accurate to second order in the time step. Although the evaluation of this term 
requires a matrix inversion, it is not demanding, as it can be performed independently for 
every block matrix corresponding to a definite angular momentum L. In addition, it is 
necessary to perform this step only once before calculating the time propagation, as this 
term is time independent. 

Calculation of the second exponential in Eq. ( ]A.3I) is more difficult, as it explicitly depends 
on time and is not block diagonal. Since and do not commute (see the structure in 
Fig. [1]), we can write this term using the split-operator formula (]A.1|) as follows: 



exp [-idlFi{t + to, At) - zdf Fx(t, At)] = exp 



--diFjit + tD,At) 



exp [-zdfFx(t, At)] exp 



--diFi{t + tD,At) 



+ O [(At)'] 



(A.5) 



The most straightforward way to calculate these exponentials is using the eigenrepresenta- 
tions of di and d5: 



exp [-idlFi{t + to, At) - idf Fx(t, At)] = U/diag jexp 
Uxdiag {exp [-zAf Fx(t, At)] } U^^U/diag |exp 



-^F,(t + t^,At) 
--^Fi{t + tD,At) 



U 



Up^.6) 
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493 where A^' are the eigenvalues of d^'^ and U/,x are the corresponding matrices of the column 

494 eigenvectors. Substitution of this result together with Eq. ( \AA\i into Eq. (1A.3P and denoting 



495 



yields 



Cj{t + At) = Vidiag < exp 



tT 
2 



diag < exp 



V2diag {exp [-zAf Fx(t, At)] } 

Vfc,(t). (A.8) 



-^F,(t + tD,At) 



496 This propagation scheme in one time step means that the wave function is first expressed in 

497 the eigenrepresentation of and propagation for a quarter of the time step is performed us- 

498 ing the field-free Hamiltonian augmented by the CAP. Then the wave function is propagated 

499 in the IR field for half of the time step. After that the wave function is transformed into the 

500 eigenrepresentation of the XUV dipole operator and the propagation for one time step is 

501 performed in the XUV field only. The wave function is transformed back to the IR eigenrep- 

502 resentation again and the propagation steps in the IR field only as well as the step using the 

503 field-free Hamiltonian are performed. This propagation scheme involves only two matrices 

504 which are time dependent and they are diagonal in different representations. It means that 

505 only two changes of representation are necessary in every time step. Note that application 

506 of Eq. flA.SP in two subsequent time steps allows for multiplication by the matrix Vf Vi 

507 between the steps. Therefore, every propagation step involves three matrix-vector multipli- 

508 cations and two diagonal matrix-vector multiplications, where only the diagonal matrices 

509 are time dependent. The matrix V2 is real and the matrix VfVi is complex symmetric. 

510 After the transformation of the initial state vector c(to) from the energy eigenrepresentation 

511 of H to the diagonal representation of d^, the rest of the propagation is performed in this 

512 basis. This approach is permissible in the context of the problem studied here, because we 

513 only need to determine the depletion of the field-free ground state at the end of the XUV 

514 laser pulse only, rather than information about the evolution of the complete wave function 
516 during the pulse. Note that only the function Fx {t, At) depends on the XUV photon energy 

516 hux and the information about the time delay to appears only in the diagonal matrix in 

517 Eq. (]A.8|) . Therefore, the same matrices Vf Vi and V2 can be used in all the time dependent 



513 calculations for all the values of hw and t^) of interest. 
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519 One more note related to Eq. (lA.SP should be made at this point. The split operator 

520 method is used there to calculate the exponential involving matrices and which do 

521 not commute. The non-zero commutator is a consequence of the approximation introduced 

522 in Sec. [TTl where we neglect the effect of the IR laser field on the ground state and neglect the 

523 XUV couplings that do not involve the ground state (see Fig. [T]). Validity of this approxima- 

524 tion implies that the commutator [d^, d^] is sufficiently small to write the exponential on 

525 the left hand side of Eq. flA.5l) simply as a product of two exponentials. This would change 

526 the form of the final time propagator (lA.SP in such way that it will contain only two time- 

527 dependent matrices. Therefore, this approach would save one matrix-vector multiplication 

528 in every time step. However, its impact on the accuracy of the time propagation needs to 

529 be tested and its implementation will be the subject of forthcoming research. 

530 In conclusion, the time required to calculate the matrices VfVi and V2 (matrix diag- 

531 onalization, inversion and the matrix-matrix multiplication) is negligible compared to the 

532 calculations of the time propagation using Eq. flA.Sp . This propagation scheme is of the 

533 second order in the time step and scales as A^^ with the dimension of the basis set used 

534 (since it involves only matrix- vector multiplications in every time step). 
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